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The phase diagram of the p-spin-interacting spin glass model in a transverse field is investi- 
gated in the limit p — » oo under the presence of ferromagnetic bias. Using the replica method 
and the static approximation, we show that the phase diagram consists of four phases: Quan- 
tum paramagnetic, classical paramagnetic, ferromagnetic, and spin-glass phases. We also 
show that the static approximation is valid in the ferromagnetic phase in the limit p — > oo 
by using the large-p expansion. Since the same approximation is already known to be valid 
in other phases, we conclude that the obtained phase diagram is exact. 
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1. Introduction 

A spin glass is a typical complex system characterized by quenched disorder and frus- 
tration. Numerous studies of spin glass systems clarified various interesting properties of 
disordered systems. 1,2 Incorporation of quantum effects into disordered systems has been of 
particular interest and studied intensively. 3-7 The noncommutativity of the operators makes 
the problem interesting but difficult, and a special care is required to obtain the correct re- 
sult. The Trotter decomposition is known to be a useful approach to treat such effects. 8 In 
this approach, order parameters become dependent on the Trotter indices and are determined 
self-consistently. 

Bray and Moore 3 proposed an approximate method to solve the problem. Their method, 
which is referred to as the static approximation (SA), is to neglect the time (or Trotter number) 
dependence of the order parameters. Using the SA, Thirumalai et al. 5 studied the Sherrington- 
Kirkpatrick (SK) model 9 in a transverse field to obtain the phase diagram. They pointed out 
the limitation of the SA, showing that the entropy does not vanish at zero temperature. 

In order to understand the nature of quantum spin glasses, it would be helpful to inves- 
tigate exactly solvable models. The infinite-range spin glass model with p-spin interactions is 
one of such tractable models. In the limit p — ► oo, this model is equivalent to the so-called 
Derrida's random energy model, 10 which can be exactly solvable in a simple way but retains 
nontrivial properties caused by quenched disorder. Goldschmidt 7 investigated this model in a 
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transverse field and obtained the phase diagram, which consists of a spin glass (SG) phase and 
two paramagnetic phases: One is the classical paramagnetic (CP) phase in which quantum 
fluctuations are irrelevant and the other is the quantum paramagnetic (QP) phase. He also 
suggested that the SA is exact in this model although there is no rigorous proof. Dobrosavl- 
jevic and Thirumalai 11 examined the validity of the SA in the same model by performing a 
systematic large-p expansion. While they showed that the phase diagram of Goldschmidt is 
correct, they also found that, for large but finite p, the SA is not valid in the QP phase. 

This model is an extreme simplification of SG models but has a great advantage to be 
exactly solvable, which has enabled us to obtain many insights not only about SG properties 12 
but also about the replica method itself 10 and information processing problems. 13 Our main 
aim in this paper is to investigate what happens in this model in a transverse field under the 
presence of ferromagnetic bias. The influence of quantum fluctuations in the ferromagnetic 
(F) phase is nontrivial, and it should be an interesting problem how the system behaves as a 
result of three competing effects of disorder, quantum fluctuations and ferromagnetic bias. In 
§2, we calculate the free energy of the model and give its phase diagram by using the replica 
method and the SA. In §3, we show the validity of the SA in the p — > oo limit by using 
systematic large-p expansion in the F phase. The last section is devoted to conclusion. 

2. Replica analysis with the static approximation 

2.1 Replica symmetric free energy 

We consider the p-spin-interacting spin glass in a transverse field. Evaluation of the par- 
tition function can be carried out by a straightforward generalization of existing methods 5-7 
to the case with ferromagnetic bias. The system is described by the following Hamiltonian: 

N 

h = - Ji,-i P <---<~rJ2< = T + v ' (!) 

ii<...<i p i=l 

where i is the site index, a z and a x are Pauli spin operators, T denotes the strength of the 
transverse field and Ji lm ..i is the quenched random interactions whose distribution function is 
given by 

jv....*> = -^)'}- w 

The limit p — > oo is taken after all calculations. The partition function of this quantum system 
can be represented in terms of a corresponding classical spin system with the Ising variable 
a = (±1) by the Trotter decomposition 8 

Z= lim Tx L-PT/m^v/m\ m = l im Zm , (3) 

where 

/ M M \ 

Z M = C M7V Trexp £ £ J^a^ . . . a ip , t + B ^ £ <7 iit a M+1 , (4) 

\ t=l i\<...<i p t=l i J 
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and the constants B and C are related to the transverse field Y as 



B = il„co.h<£, C-fi-tah^V. (5) 
2 M \2 M J y J 

The symbol Tr denotes the trace over the cr-spins. 

We use the replica method: 9 

\Z n ] — 1 

[log Z) = lim [ - — i (6) 

ra-»0 n 

to carry out the random interactions Ji x ...% v averages [•••]. The replicated partition function 
is given by 

[Z M ] = Trexp hj^T^ ^ hvZ.<^' +ifLLbMt 

\ t,t'=l/J.,v=l V i / t=l /i=l \ i 

M n 

+ 5 EEE<<m)> (7) 

t=l /i=l i 

where the replica indices are denoted by \x and v. We have omitted some irrelevant constants. 
The spin product term (^2i a it a it' can ' 3e simplified by introducing an order parameter 
q^ and its conjugate Lagrange multiplier for the constraint q^ = J2i a it a it'/^ '• ^ n ^ ne 
present case, we must distinguish diagonal q^ and off-diagonal q^(n / v). Physically, q^jf is 
a measure of quantum fluctuations and q^, is the SG order parameter. If there are no quantum 
fluctuations, the spin configuration af t does not depend on time t and q^ = 1. Quantum 
fluctuations reduce from unity. Hence, 1 — q+lf gives a measure of quantum fluctuations. 
Also, to simplify ^X^^ft/^) > ^ ne ferromagnetic order parameter is introduced and the 
constraint = ^2i&i t /N is imposed by the integration over the conjugate variable m^. 
Using these notations, we can rewrite the effective partition function as 

™ = fun d < d ^ n n n n 

x exP J£ £ (££ «r - ^w) + EE W)' - i 5 "'*') 

+ EE(t W)P "¥™ W ) +logTrexp(-^ cff )|, (8) 

where 

^ = - fi EE ^^+1 - ^ E E -^EE <' - jp E E 

t A* M t t,t' M t^t' 

(9) 

We can calculate the free energy per spin F of the replicated systems in the thermodynamic 
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limit by the saddle-point method. The result is 

= E E (IS ^-) P " ) «) p - 

t,t' H<v X 7 t,t' ft V 7 

+ E E (§r K) p - + lo § Tr ex p ' ( 10 ) 



where 



g# = <«) , = \pJ 2 P { q %Y-\ (lib) 



r = <of> , m^ = (3j 0P (m^r-\ (lie) 



The brackets (• • • ) denote the average by the weight exp(— H e g). It is difficult to solve these 
equations exactly for arbitrary values of p because of the time dependence of the order param- 
eters. However, in the limit p — ► oo, the problem is considerably simplified because conjugate 
variables can be either or oo. For example, if we restrict ourselves to the case that m^ is 
non- negative (the other case can be treated similarly), eq. (11c) implies < < 1, which 
leads to either — ► (0 < < 1) or — > ±oo {m^ — > 1). Therefore, the SA gives the 
exact phase boundaries as is shown below. 

Now, we assume the replica symmetry (RS) and use the SA: 

(ftv = R ' 1w = 1 ' m t = m - ( 12 ) 
As already noted, R is the order parameter measuring the effect of quantum fluctuations, q is 
the conventional SG order parameter and m denotes the ferromagnetic order parameter. The 
free energy is then reduced to the expression 

—(3f = — lim lim — — 
M— >oon— >0 n 

= -J 2 J 2 (R p - q p ) + -qq - RR + fij m p - mm + lim lim - log Tr exp (-H eS ) . (13) 

Since all order parameters and the conjugate variables are independent of time and replica 
indices, the summation of spin products in H e g is rewritten as 

E E - \ \ { E E «r) 2 - E (E «f) 2 ) • <«) 

n+v t,t' ( \ (i t / n \ t / J 

Then, the effective Boltzmann factor exp(— H c q) is considerably simplified by the Hubbard- 
Stratonovich transformation. We obtain 

exp(-# cff ) = j D Zl J] |y Dz 2 exp (b £ a?a? +1 + A £ ^ | , (15) 

where 



A= J2R-qz 2 + y/qzi + m , Dz=-^Le~5 22 . (16) 

V 2vr 
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Using the Trotter formula, we can take the spin trace in the limit M — > oo as 

lim Tr expfflV a t a t+1 + ^- Y] <*t | = Tie AaZ+l5TaX = 2 cosh \J A 2 + /3 2 r 2 . (17) 

V t t / 

Because each replica gives the same contribution to the replicated free energy, the limit of 
n — ► is easily taken. The result is 

-Pf = i/3 2 J 2 (R p - q p ) + hqq -RR + (3j m p - mm + J Dz 1 log J Dz 2 (2 cosh u), (18) 



where 

|2 , o2-n2 



w = (A + (3 T Y . (19) 
The saddle-point conditions of the free energy are 

m = pjopmT 1 , (20) 

R = ^(3 2 J 2 pRP-\ (21) 

q= l -(3 2 J 2 Pq v-\ (22) 

m = J Dz x Y~ x J Dz 2 Auj- x sinh u, (23) 

i? = y Dz 1 Y~ 1 [J Dz 2 A 2 uj- 2 cosh uj + /? 2 r 2 y L>z 2 ^~ 3 sinhL^ , (24) 



g = y Dz x Y~ 2 (^J Dz 2 Au' 1 sinh u?j , 



(25) 



where 

y = j Dz 2 coshw. (26) 

When r is equal to 0, R becomes unity and the free energy is reduced to the classical result. 10 ' 12 
These eqs. (20)- (25) generalize the result of Goldschmidt 7 to the case of finite jo- 
2.2 Solutions of the equations of state. 

The following inequality is useful to evaluate the solutions of eqs. (20)-(25). 14 

R = j Dz x Y- x (J Dz 2 A 2 u>~ 2 cosh u + f3 2 T 2 j Dz 2 cu' 3 sinh u?j 

> y Dz 1 Y~ 1 y Dz 2 A 2 lj~ 2 cosh ca> J Dz 1 Y~ 1 J Dz 2 A 2 aj~ 2 sinhu; 

> J D Zl Y- 2 (^j Dz 2 Auj- 1 sinh u^j =q. (27) 

From eqs. (20)-(22) and < R, q, m < 1, conjugate variables can be either or oo in the limit 
p — ► oo. These conditions restrict combinations of the values of conjugate variables to 

(fh, R,q) = (0, 0, 0), (0, oo, 0), (0, oo, oo), (oo, 0, 0), (oo, oo, 0), (oo, oo, oo). (28) 
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We must exclude the solution (fh,R,q) = (oo,0,0) as follows. Substituting these values 
(oo,0,0) into eqs. (23)-(25), we find the solution (m,R,q) = (1,1,1), which contradicts the 
condition (R, q) = (0, 0) as one can check from eqs. (21) and (22) with p — > 00. Similarly, 
(fh,R,q) = (oo,oo,0) is also inappropriate. Finally, we get four solutions 

(m, R, q) = (0, tanh/3r, 0), (0, 1, 0), (0, 1, 1), (1, 1, 1). (29) 

The solution (m,R,q) = (0, 1,0) is a paramagnetic one which is identical to that in the 
case without quantum fluctuations. In this solution quantum effects are irrelevant so that this 
phase is the CP phase. The solution (m, R, q) = (0, tanh (3T, 0) represents a non-trivial 
paramagnetic phase which arises due to quantum fluctuations. R is reduced from unity by 
the transverse field T and this phase is the QP phase. The solution (m, R, q) = (0, 1, 1) is for 
the SG phase. In the limit p — > 00, the finite value solution of q is limited to q = 1 from eqs. 
(20)-(22). From the inequality R > q, R is also equal to 1 and hence we expect that quantum 
fluctuations are irrelevant in this phase. 11 ' 15 The solution (m, R, q) = (1, 1, 1) is ferromagnetic. 
In this phase, all order parameters are restricted to 1 in the limit p — > 00. Therefore, quantum 
fluctuations are also irrelevant in this phase as will also be shown below. We summarize the 
above results in Table I. 



Phase 


(m,R,q) 


(m,R,q) 


CP 


(0,1,0) 


(0,oo,0) 


QP 


(0,tanh/?r/(/?r),0) 


(0,0,0) 


F 


(1,1,1) 


(00, 00, 00) 


SG 


(0,1,1) 


(0, 00, 00) 



Table I. Values of order parameters in various phases. 



Substituting the values of the order parameters and conjugate variables into eq. (18), we 
obtain the corresponding free energy as 

/ QP = -T log 2 - T log cosh (3T, (30) 

/cP = -^</ 2 -Tlog2, (31) 

/f = -jo, (32) 

fsG - -tJ— - -00 (p - 00). (33) 

V 7T 

For this large p case the free energy of the SG phase is always smaller than the other free 
energies in eqs. (30), (31) and (32), although one knows that in spin glass phases maximization 
is the correct procedure. 1 However, this feature is an artifact of the RS solution. 
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2.3 Free energy in the spin glass phase 

The correct solution of the SG phase is derived by the RSB. In the limit p — ► oo, the first 
step of the RSB (1RSB) is known to be sufficient in the classical case (T = 0). 12 We may 
expect that the same is true in the presence of a transverse field and therefore discuss the 
1RSB scheme here. Note that at sufficiently low temperatures the full-step RSB is known to 
be necessary for finite j^ 16 ' 17 which we do not take into account explicitly here because we 
take the limit p — > oo in the end. New order parameters and a branch-point parameter mi 
are defined as follows 

q ¥iM = qo ; = Ql \^\^ (34) 

where I = 1, . . . , njm\ is the block number and fj,i = 1, . . . , m\ is the index inside a block. 
Detailed calculations are given in Appendix A. The 1RSB free energy is found to be given as 

Pf = Jmi ( b 2 J 2 (qo) p - qoqo) + U 1 ~ ™i) (^AqiT ~ QiQi 



-^(3 2 J 2 (R) P + RR- (3j m p + mm-^- j Dzi log J Dz 2^j Dz 3 2 cosh u^j (35) 

where 

wi = {A\+(3 2 T 2 Y , A x = y/% Zl + - q z 2 + \J 2R - q\ z 3 + m. (36) 

If we set qo = qi, the RS solution (18) is reproduced. The saddle-point conditions of the 1RSB 
free energy are given by 

m = (3j pm p -\ (37) 

R = ^(3 2 J 2 pRP-\ (38) 

90 = \(i 2 J 2 vq V o-\ (39) 
qi = \p 2 J 2 pq P ~\ (40) 
m = J DzxYf 1 J DZ2Y™ 1 - 1 , J Dza^iu^sinhwi, (41) 

R = j Dz x Y x - x J DZ2Y™ 1 - 1 (^j Dz 3 A\uj x 2 cosh wi + (5 2 T 2 j Dz 3 u;7/ 3 sinh u^j ,(42) 
q = f Dz i (y' 1 J DZ2Y™ 1 - 1 J DzsA^smhoj^J , (43) 
qi = J Dz! yf 1 J Dz 2 Y 2 mi_2 Dz 3 ^l^ 1 sinhc^ , (44) 

where 

Yi = Fi(zi) = y Dz 2 ^ £>*3 cosh^ , Y 2 = Y 2 (z 1 ,z 2 ) = j Dz 3 coshwi. (45) 

Inequalities i? > gi > qo are also derived in a similar way to the derivation of eq. (27) 
as explained in detail in Appendix A. In the limit p — > 00, we find from these inequalities 
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and eqs. (37)-(44) that the only non-trivial RSB solution is (m,R,qo,qi) = (0,1,0,1). The 
corresponding free energy is 



f = --BJ 2 mi -T— log 2. 
4 mi 



Taking a variation with respect to mi, we find 



/sg = -J\/\og2. 



(46) 
(47) 



This is the correct free energy of the SG phase. 

We can determine all the phase boundaries by equating free energies between different 
phases. The phase diagram thus obtained is shown in Figs. 1-4. We can see that, as V grows, 
quantum fluctuations reduce the ferromagnetic order and cause a phase transition to the QP 
phase as in Fig. 2. Order parameters discontinuously change at any phase boundary. In that 
sense, all the phase transitions are of first order. 




Fig. 1. Full phase diagram of the model in the 
limit p — ► oo. The QP phase appears when 
r/J becomes larger than Ti/J — l/\/2. The 
CP phase is completely suppressed by quan- 
tum fluctuations in the range ofr/J>r f /J = 
log(2+v / 3) / (2-^/2) ■ Replica-symmetry breaking 
exists only in the SG phase. 



r/J 



1/2 

(In2) 



(In2j 



/2 



T=0 









QP 








F 




SG 







Jn/J 



Fig. 2. Ground-state phase diagram on the T = 
plane of Fig. 1. For large T, quantum fluc- 
tuations destroy the ferromagnetic order and 
cause a phase transition to the QP phase. 



3. Validity of the static ansatz 

3.1 Expansion from the large-p limit 

In this section, we check the validity of the SA. Again, the method is generalization of 
ref. 11 to the case with ferromagnetic bias. For that purpose, we introduce corrections to the 
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QP 






CP 






SG 


F 





(In2) 



1/2 



Jn/J 



T/J 




Fig. 3. Schematic phase diagram on the T-jo 
plane in the range of Ti < T < Tt- For large 
T, the QP phase appears instead of the CP 
phase. As T grows, the CP phase diminishes 
and vanishes at r = Tt- 



Fig. 4. Phase diagram on the T-jo plane for large 
T. The SG phase disappears at T/J — T /j = 



SA with t, i'-dependence and expand the free energy with respect to those correction terms. 
Then, it is shown that the time-dependent parts are irrelevant in the limit p — > oo. 
We start from the RS solution 

g$ = R(t, t') , qj£ = q(t, , < = ro(t). (48) 

Separating each conjugate variable to the static and time-dependent parts, we rewrite the 
effective Hamiltonian H e R = H sta t + V(t,t') as 

*«« = -b E £ « +1 -^EE ^ - ^E E « - E E ^ (49) 
ffr o = -^EE A ^ *W -^EE *'K%< - ^ E E A ^K< ( 5 °) 

M t^t' ni=v t,t' M * 

where 

R(t, t')=R + AR(t, t') , q(t, t') = q + Aq(t, t') , m(t, t') = rh + Am(t, t'). (51) 

It is expected that the order parameters are monotone decreasing functions of the time interval 
\t — 1'\ because they are originally written as spin-correlation functions eqs. (lla)-(llc). Then, 
the time-dependent parts of their conjugate variables, which are the pth powers of the order 
parameters, become drastically small for large p. Therefore, it is reasonable to expand the free 
energy with respect to the time-dependent part. The free energy is expanded with respect to 
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V(t, t') to first order as 

M = \ W E (V A(t, t'y - - A ^> t')) 

t,t> v ' 

-TF E (Plom(t) p ~ mm(t) - Am(f)m(t)) - lim 1 (logTrexp(-tf stat ) + (V% at ) . (52) 

The brackets (• • • ) s t a t denote the average by the weight exp(— H sta t). The equation satisfied 
by each order parameter can be obtained by taking a functional derivative with respect to the 
time-dependent part of the conjugate variable. The results are 

™(*) = ^\E^) > ( 53a ) 

\ P / stat 

q(t,l/) = -\imUj2^t) , (53b) 



n— >0 n 



stat 



^,0 = lhni(E«) • ( 53c ) 

\ ^ / stat 

The new Boltzmann factor exp(— -ff s t a t) is identical to exp(— H e g) under the SA and can be 
simplified by the Hubbard-Stratonovich transformation as eq. (15). From eqs. (15) and (53a)- 
(53c), we can verify that q and m are time-independent and only R is time-dependent because 
of the independence of each replica in i? s tat and the translational invariance in the Trotter 
direction. Then, the problem is reduced to the evaluation of the correlation function of the one- 
dimensional Ising system in a field. Calculations are somewhat involved but straightforward. 
Details are given in Appendix B. The result is 



(54) 



R(t) = J Dz x Y- 1 (^j Dz 2 A 2 uj' 2 cosh u + (3 2 T 2 J Dz 2 uj~ 2 coshcj(l - 2t)\ . 

We have used the same notation as in eqs. (16)- (19) and the continuous-time notation 

t - t' 

t = lim . (55) 

M->oo M 

The result (24) of the SA is reproduced by integrating eq. (54) over r. The equations of other 
order parameters are identical to the result under the SA. Hence, from eq. (27), we see that 
the inequality R(t, t') > q still holds for any r. 
3.2 Explicit solution of R(r) in the ferromagnetic phase 

In the F phase, all conjugate variables go to infinity in the limit p — > oo and it is reasonable 
to assume 2R = q according to eqs. (21) and (22). Hence, the integration over z 2 just gives 1 
and we find 

R F (r) = [DzJ A 2 ^ 2 + ^ r2tJ -2 CosMl-2T) \ 
/ V coshu; / 
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The conjugate variables fh and q are proportional to p and very large which enables us to derive 
the leading finite-p correction by the systematic large-p expansion. 11 ' 15 For large conjugate 
parameters, we can estimate the integral (56) by the saddle-point method for fixed r. To 
compare the time-dependent result with that under the SA, we start from the result of the 
SA. Under the condition 2R = q, eq. (24) reads 

R = j D Zl (a 2 uj- 2 + (i 2 T 2 J Dz 2 " _3 |^) ■ (57) 

The first term on the right-hand side of eq. (57) is rewritten as 

f oof (l + y/qzjm) 

/ Dzi A 2 oj- 2 = / Dz x ^ '— (58) 

J J l + 2^fqz l /m+(^/qz l /m^ +{(3T/m) 2 . 

For large p, because 1/rh oc 1/p is very small, we can expand the right-hand side of eq. (58). 
After straightforward calculations, we obtain the leading 1/p-correction term as 

D Zl A 2 u- 2 « 1 - 5f = 1 - (- V \m- 2 ^ 2 . (59) 



/ 



m2 \hj P 2 



Similarly, the second term on the right-hand side of eq. (57) is rewritten as 



c 



-2ui 



D Zl (^/3 2 r 2 u;- 3 - 2/? 2 r 2 o;- a 1 ^ e _ 2a) j , (60) 

and exp(— 2uj) is exponentially small for large p. The first term of eq. (60) can be evaluated 
by the series expansion with respect to 1/p as in eq. (59). The result is 

This term is at most proportional to 1/p 3 and is negligible. We can also evaluate other order 
parameters m and q. Under the condition 2R = q, eq. (23) is rewritten as 

m = J Dz\Y~ x J Dz2 Auj^ 1 sinhw = j Dz\ Aj^tanhu; 



/ 



D Z1 Au- 1 « 1 - \ (-) 2 - 2 m- 2 ^ 2 « 1 - \ ( IV 1. (62) 
2 \JoJ P 2 2 \j J p 2 



Similarly, eq. (25) reads 

q = J Dz x Y- 2 (^J Dz 2 Auj- x sinhu;^ = J Dzi (Ajtanhu) 2 

« j D Zl A 2 u- 2 ^l-{^j ^ 2 =R. (63) 

This equation is consistent with the condition 2R = q. 

Next, we proceed to a time-dependent analysis. From eq. (56), its first term on the right- 
hand side gives the same corrections as in the SA case. We rewrite the second term as 

/" 2 coshw(l-2r) /" 2 e~ 2 ™ + e"^ 2 " 2 ^ /" 2/2™ ^22^ 

J 1 U coL = J 1 TTe- 2 ^ * J DZ ^ ( 6 +^ (2 ' 2T) ) • 

(64) 
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This integration can be evaluated by the saddle-point method. The saddle-point condition of 
oj~ 2 exp(— 2tlo) is 

, 2 = 4r^fl-^V (65) 



It is difficult to solve this equation exactly. However, for large p, the second term on the 
right-hand side of eq. (65) is vanishingly small and we can approximate the saddle point as 
z = ±2r-\/(^ The contribution from the saddle point z = 2ry / g is 

|(2Tg + m) 2 + /3 2 r 2 } _1 exp ^-2r 2 g - 2t\J {2rq + mf + p 2 T 2 ^j 



" ^2 J2 «P |- [r 2 J 2 (3 2 + 2r^ (r^^ + (3j ) + ^ j pj . (66) 

The other saddle point z = —2T^fq gives a similar contribution. Replacing 2r by 2 — 2r, we 
can also obtain the saddle-point value of w~ 2 exp(— uj(2 — 2r)). Then, the explicit result of 
Rir) in the F phase is given by 

R F (r) « 1 - f-V ^ + f-1 ^/(r,p) (67) 
VW P 2 VW P 2 

where the function f(r,p) expresses the time-dependent correction which decreases exponen- 
tially as p grows. Because the third term of eq. (67) is vanishingly small, the main finite-p 
correction is the second term, which is identical to the SA result. Accordingly, in the F phase, 
the time-dependent part of the finite-p correction is exponentially small for large p as in the 
CP and SG phases 11,15 and the SA is valid in that sense. We also calculated the free energy to 
the order I /p. However, we found that the first order correction vanishes and the free energy 
remains as in eq. (32). To obtain the leading order in 1/p, we must proceed to the next order 
approximation but it is beyond our purpose in this paper. 
3.3 Remarks 

In the F phase, we have found that the leading corrections of order parameters are actually 
time independent. For other phases, previous works revealed that in the CP and SG phases 
similar results hold. 11 ' 15 Hence, in these instances, the SA is valid not only in the limit p — ► oo 
but also as long as p is adequately large. Meanwhile, in the QP phase, strong disagreement 
occurs between the SA and the time-dependent analysis for large but finite p. 11 Not only the 
spin autocorrelation function R(t) is actually time dependent, but also the low temperature 
behaviour shows violations of the thermodynamic law within the SA. From these facts, we 
may conclude that the SA well describes the region in which quantum fluctuations are weak, 
but strong quantum effects lead to a collapse of this approximation and unphysical behaviours 
of thermodynamic quantities. In spite of such inexpediency, the free energy recovers the SA 
results in the limit p — ► oo even in the QP phase. Consequently, the SA gives correct free 
energies in all the phases in the limit p — > oo, and the phase diagram depicted in Figs. 1-4 
should be exact. 
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4. Conclusion 

In this paper, we have studied the p-spin-interacting spin glass model with ferromagnetic 
bias in a transverse field by the replica method. Trotter decomposition has been employed to 
reduce the quantum system to a classical one and the SA has been assumed to obtain the 
solutions of the equations of state. We have clarified the structure of the full phase diagram, 
which consists of four phases, the CP, QP, F, and SG phases. 

We have also checked the validity of the SA in the F phase by the large-p expansion. 
Leading finite-p corrections of the order parameters have been calculated and it has been 
shown that they are actually time-independent. It is known that similar results hold in the 
CP and SG phases. This is not the case for the QP phase. Nevertheless the free energy of the 
QP phase turns out to be identical to the SA result in the limit p — > oo. In that sense, the 
SA gives correct solutions in this limit and our phase diagram is exact. 

Admittedly, the model investigated in this paper is not a faithful reproduction of real SG 
systems. However, it has a great advantage that order parameters and the free energy can be 
exactly obtained. Goldschmidt found two qualitatively different types of paramagnetic phases, 
CP and QP. We also found that quantum fluctuations reduce the ferromagnetic order and 
cause a transition to the QP phase. These properties appear to be plausible in more realistic 
systems. On the other hand, we saw that in the CP, F, and SG phases quantum fluctuations 
are completely irrelevant, which should be specific to this model. For more realistic models 
(like the SK model) we should take into account the influence of quantum fluctuations on 
order parameters. However, it is difficult to treat such an effect as was done in the present 
paper and we need different techniques. Finding effective approaches and improving the SA 
remain interesting problems to be investigated in the future. Our present results will serve as 
a first step to understanding the interplay between quantum fluctuations, ferromagnetic bias 
and quenched disorder. 

After submission of the manuscript, we came to notice that the same problem was discussed 
by Saakyan 18 and Jun-Ichi Inoue. 19 The originality of our work lies in the systematic analysis 
of the validity of the SA and the explicit clarification of the structure of the full phase diagram. 
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Appendix A: 1RSB free energy and the inequalities R > qi > qo 

In this Appendix we derive the 1RSB free energy eq. (35) and the inequalities R > q\ > qo- 
Under the 1RSB scheme and the SA, eq. (10) reads 

Pf = \™i Q/3 2 J\l - §b9o) - \(mi - 1) Q/3 2 J 2 q p - gift) - (^f3 2 J 2 R P - RR^j 

— (Pjom p — mm) — lim lim — log Tr exp (— H c g) • (AT) 

M-+oo n^O 71 

The effective Hamiltonian H c g is written in the form 

- HcS = b £ e +™EE^ +i *EE « 

/it M i /U t^t' 

/ n/mi block \ 

+1 » E E + (£ -a) E E E « • ( A - 2 ) 

V l Vl+vi t,t' j 

We can rewrite the summation of spin products as 

-h cS = \m ( E E 2 + \^~ ft) E ( f E E °?) 2 - E f E < 

\ n t / i [ V m t / \ i 

+ (* - E f E <*) + 5 E E + - E E (a-3) 

To take the spin trace, the Hubbard-Stratonovich transformation is employed for the quadratic 
terms. The result is 

n/mi ( lm\ 



e -H e g 



= J Dz x H I J Dz 2 H (J Dz 3 e L ) 
l «;=(7-l)mi+l 



(A4) 

m=(l-l)mi+l 

where L = Ax/MJ^t^t' + B J2t a t l a tU and A i is defined in ec l- ( 36 )- Using the Trotter 
formula, we can take the limit M — > oo and perform the spin trace as in eq. (17). The result 
in the limit n — > is given as 

lim lim - log Tr exp (-H e g) = — [ Dzi log f Dz 2 ( [ Dz 3 2coshwi N ) . (A-5) 
M^oomOn mi J J \J J 

The equations of state eqs. (37)- (44) are obtained by taking a variation of the free energy eq. 

(35) with respect to conjugate variables and order parameters. 

Next, we derive the inequalities R> qi > qo- From eqs. (41)-(44), 

R = j Dz± Ff 1 J Dz 2 Y™ 1 ' 1 (^j Dz 3 A\u^ coshwi + (3 2 T 2 J Dz 3 u^ 3 sinhwi) 

> J Dzi Ff 1 j Dz 2 Y™ 1 ' 1 J Dz 3 A\u^ 2 coshc^ 

> j Dzi Yf 1 J Dz 2 Y~ 2 mi_1 J Dz 3 A\u^ 2 sinh^i 

> j DziYf 1 J Dz 2 Y^ 2 (^j Dz 3 A l u^ 1 smhu^ = q x . (A-6) 
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Similarly, we can show q\ > go from the definition of Y\ (45) as 

gi = J Dz 1 Y 1 - 1 j Dz 2 Y^ 1 - 2 (^J Dz 3 A^ 1 smhu^ 

> J D Zl ^Yf 1 J Dz 2 Y™ 1 - 1 J Dz 3 A^ 1 sinhu^ = q . (A- 7) 

Appendix B: Evaluations of the correlation function 

We calculate the correlation function and derive the expression for R(t) eq. (54). The un- 
normalized correlation function of the one-dimensional Ising system with periodic boundary 
is 

(M M \ 

jJ^Wt+i + hJ^vt • (B-l) 
t=i t=i J 

We can compute the correlation function G(t, t') by the transfer matrix method. The general 
solution is 

G(t, t') = 4x 2 + x 2 _ [X 1 - 1 ' X M ~ {t ~ t,] + \+ - (t - t ' ) A t _- t ') + (2x 2 + ~ + ( 2x - ~ ( B " 2 ) 

where A± are the eigenvalues of the transfer matrix and x± are the first components of the 
eigenvectors |±). Their explicit forms are 

A ± = e J (cosh h ± \/cosh 2 h-l + e" 4J ) , (B-3) 

!±) = D± ( e J {smhhT Vsinh 2 + e- 4J ) ^ = ( y ± ) ' (B4) 
where Z)± are normalization constants. In the present case, the parameters J = B = 

1 /2 

log (coth/?r/M) ' and h = A/M depend on the number of spins M. Restoring the omit- 

1 /2 

ted overall factor C = {(1/2) sinh2/3r/M} ' , we can obtain the finite value of the correlation 
function in the limit M — ► oo. After straightforward calculations, we get 

where u> is given in eq. (19). Hence, the correlation function G(t,t') is given by 

G(t, t') -► G(t) = 2^ 2 lj- 2 coshw + 2(/?r) 2 Lj- 2 coshw(l - 2r). (B-6) 
Substituting eq. (B-6) into eq. (53c) and taking the limit n — ► 0, we finally get 

R(t) = J Dz ± Y- 1 (^J Dz 2 A 2 oj- 2 cosh u + (3 2 T 2 J Dz 2 to' 2 coshu;(l - 2r)^ . 
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